library(forecastLM)
md1 <- trainLM(input = ny_elec,
               y = "y",
               seasonal = c("month"),
               trend = list(linear = TRUE))


plot_res(md1)
md2 <- trainLM(input = ny_elec,
               y = "y",
               seasonal = c("month", "wday"),
               trend = list(linear = TRUE))


plot_res(md2)
md3 <- trainLM(input = ny_elec,
               y = "y",
               seasonal = c("month", "wday", "hour"),
               trend = list(linear = TRUE))


plot_res(md3)
md4 <- trainLM(input = ny_elec,
               y = "y",
               seasonal = c("month", "wday", "hour"),
               trend = list(linear = TRUE),
               lags = c(1:24))


plot_res(md4)
md5 <- trainLM(input = ny_elec,
               y = "y",
               seasonal = c("month", "wday", "hour"),
               trend = list(linear = TRUE),
               lags = c(1:24, 48, 72))


plot_res(md5)